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Abstract 

Despite attention focused on tlie population status and rebuilding trajectory of Atlantic bluefin tuna {Thunnus thynnus), the 
reproduction and spawning biology rennains poorly understood, especially in the NW Atlantic. At present, the eastern and 
western spawning populations are believed to exhibit different reproductive characteristics and, consequently, stock 
productivity. However, our study suggests that the two spawning populations, the Gulf of Mexico and the IVlediterranean 
Sea, could show similar reproductive features and spawning strategies. Between 2007 and 2009, gonad samples from 
female Atlantic bluefin tuna were collected in the northern Gulf of Mexico (n = 1 47) and in the western Mediterranean Sea 
(n = 40). The histological and stereological analysis confirmed that sampled eastern and western bluefin tuna exhibit the 
same spawning duration (three months) but the spawning In the Gulf of Mexico begins one month earlier than in the 
Mediterranean Sea. Western bluefin tuna caught in the peak of the spawning season (May) showed a similar spawning 
frequency (60%) to the spawning peak observed in the Mediterranean Sea (June). Fecundity for the Gulf of Mexico fish 
(28.14 eggs g^') was lower but not significantly different than for fish sampled In the Mediterranean Sea (45.56 eggs g^'). 
Our study represents the first comparative histological analysis of the eastern and western spawning stocks whose findings, 
combined with new determinations of size/age at maturity and possible alternative spawning areas, might suggest basic life 
history attributes warrant further scientific and management attention. 
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Introduction 

The reproductive biology of Atlantic bluefin tuna (ABFT; 
Thunnus thynnus, L. 1758) remains poorly understood despite the 
high economic value of this fishery and its exploitation throughout 
the Atlantic Ocean and the Mediterranean Sea. These uncertain- 
ties directly afiect our understanding of the recruitment and 
productivity of the stock and could result in inefficient manage- 
ment of the fishery. The International Commission for the 
Conservation of Atlantic Tunas (ICC AT) manages the population 
as two stocks (eastern and western), separated at 45°W based on 
two known spawning areas, the Gulf of Mexico (GMX) and the 
Mediterranean Sea (MED). This division assumes ABFT exhibit 
natal homing [1,2] and have different maturation schedules [3,4]. 
Nevertheless, electronic tagging and genetic/chemical marker 
studies have shown stock mixing on the foraging grounds 
indicating more complex population dynamics [1,2,5-7]. The 
eastern stock is estimated to be tenfold larger than the western 
stock [8] ; consequently, the mixing rates are unbalanced with the 
eastern stock having greater influence on the western population, 
and any management action aimed at the eastern stock may 
indirectly afiect the western stock [2]. Understanding the 



reproductive potential of both stocks is essential as it influences 
recruitment and hence the sustainabUity of the stocks and their 
capacity for supporting commercial fisheries [3,9]. 

Electronic tagging and macroscopic examination of gonads are 
useful tools but, alone, are not sufficient for assessing population 
reproductive dynamics. Studies based on the histological analysis 
of gonads have allowed the characterization of the reproductive 
cycle of eastern spawners by the identification of maturation stages 
and estimations of several reproductive parameters [4,10-13]. In 
the western and central Adantic, these studies are scarce or 
impaired by large uncertainties [3,9,14]. Past studies examined 
gonad histology in fish from the Bahamas and the mid-Atlantic 
Bight and found a significant proportion of post-spawning females, 
but no ovaries with ovulated oocytes were observed [9,15,16]. 
However, these studies provided limited information since the 
samples were taken far from the known spawning areas and could 
not provide information about the reproductive condition of fish 
classified as active spawners [17], such as spawning frequency 
and/or fecundity [18]. 

While age at maturity in the eastern stock has been estimated as 
3-4 y [4, 1 0] , this parameter for the western spawning stock is the 
object of intense debate. ICCAT assumes an age at 50% maturity 
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for the western stock of 8 y, but other studies provide estimations 
ranging from 4—16 y [9,14,19-21]. The younger maturation age 
observed in the eastern population [4,9-11] leads to higher 
lifetime reproductive output and, consequendy, larger productivity 
in the eastern than in the western stock. 

Fecundity estimates allow the quantification of the reproductive 
capacity of individual fish and are essential for accurate assessment 
of the spawning stock biomass [22]. Assuming environmental 
characteristics between the MED and GMX spawning grounds 
are different, and ABFT exhibit natal homing, fecundity must be 
calculated separately for each stock. 

An intensive evaluation of the maturity status of ABFT in the 
GMX compared with the MED would improve our comprehen- 
sion of the reproductive connections between both stocks. In the 
present study, a histological and stereological analysis of gonads 
from ABFT caught in the GMX and the MED was undertaken to 
1) determine the reproductive status of female ABFT in the GMX 
through qualitative and quantitative histological studies, and 2) 
compare the reproductive parameters estimated between both 
spawning grounds. 

Materials and Methods 

Sample collection and biometry 

Female ABFT were sampled from commercial fisheries in the 
GMX and the MED between 2007 and 2009. The National 
Marine Fisheries Service Pelagic Observer Program sampled fish 
from longline vessels in the northern region of the Gulf of Mexico 
from February^uly (« = 147). Eight samples obtained in February 
(n = 1) and March {n = 7) were not included in the stereological 
analysis because of low monthly sample size. In the MED, ABFT 
were sampled from longline fishing vessels on the western MED 
spawning grounds from mid June to mid July, 2008 (m = 40). For 
all samples, curved fork length {CLp) of each individual was 
measured to the nearest cm and converted to straight fork length 
(Lp) using the formula: Lf = CLf -0.955 [23]. Body mass {Mg) 
was calculated from Lp by location and timing of catch according 
to ICCAT conversions [23]. Ovaries were immediately removed 
and weighed, and the volume of the pair of ovaries (Fo) was 
estimated from their mass {Mq) according to the equation 
Fo = 0.9 174 -Mo [12]. The gonadosomatic index (Iq) was 
calculated from /g = lOQ-A/o-A/^'. 

Histology 

A subsample (0.5-1 cm^) was removed from the central portion 
of one ovary and fixed for at least 24 h in 4% formaldehyde (lO'l^i 
formalin) in phosphate buffer, 0.1 M, pH 7.2. Tissue samples were 
dehydrated through increasing concentrations of ethyl alcohol, 
cleared, and embedded in paraffin wax. Samples were cut into 5-6 
fim sections using a microtome, and stained with either 
haematoxylin-eosin (GMX) or haematoxylin-VOF (MED) [23]. 
All images were taken on a Leica DMI 6000b using Image-Pro 
Plus. 

Based on the oocyte development, eight distinct types of 
developing follicles were distinguished: perinucleolar (PNF), lipid- 
stage (LSF), vitellogenic (VF), and oocyte maturation follicles 
(OMF), which consisted of migratory-nucleus (MNF) and hydrated 
(HF) follicles. Additionally, a- and jS-stage atretic follicles (aAF and 
)SAF, respectively) and postovulatory follicles (POF) were counted 
(Figure 1). Based on the most advanced follicle type in the ovary 
and the extent of atresia, fish were classified as inactive (IN), active 
non-spawning (ANS), or active spawning (AS) [17,25]. 



Spawning frequency 

Spawning frec]uen('\" ^\'as estimated by the posttn'ulatorv' follicle 
method [26] as adapted by [27]. This method calculat(;s the m(;an 
spawning fraction as the total number of spawning females whose 
ovaries show postovulatory foUicles (POF) divided by the total 
number of mature females sampled. 

Stereology 

A stereological model-based method was applied to estimate 
numbers of the diflFerent categories of foUicles according to the 
formula, NV = KP^^ N]f'V~'^-^ where Nv is the numerical 
density (number per unit volume) of the considered follicle typt, 
/? is a shape coefficient, A'is a size distribution coefficient, Na is the 
number of follicle transections p(;r unit area, and Vy is the partial 
area (volume density). Vy was calculated by analyzing 10 digital 
micrographs from each ovary using ImageJ [18,28]. Values of 
both P and A" previously calculated for ABFT were used for this 
study [11,12,29]. The total number of foUicles (jV) was calculated 
by extrapolating Nv to the total ovarian volume, N = N yVo- 
Ovarian volume loss through processing was measured, and 
corrections were applied, 34.8% for MED and 43.3% for GMX. 
Finally, the relative number of follicles (number of foUicles per 
gram) was estimated as Ng = N'Mj^^ . 

Statistical analysis 

Comparisons of means of the stereological and biometrical 
parameters among years were performed using the Kruskal-Wallis 
test, and parameters with no significant difference were regrouped 
by month (a = 0.05). Monthly variation was also analyzed using 
the Kruskal-WaUis test (a = 0.05). The Mann-Whimey U-test with 
Bonferroni correction was used to assess significant cUfferences 
between pairs of months (a = 0.0125) [30]. 

Results 

The Lp for fish sampled in the GMX (n=147) was 172- 
326 cm, and 120-240 cm for fish sampled in the MED (n = 40). 
The Ig was 0.32-6.9 in tiie GMX and 0.3-5.8 in tiie MED. Mean 
Lp±SD offish sampled in tiie GMX (235.61 + 19.81 cm) was 
significantiy larger (Mann-Whitney, p < 0.001) than those sampled 
in the MED (198.15 + 27.84 cm; Table 1). No significant 
differences in Lp were observed within the GMX throughout 
the sampling period (Kruskal-WaUis, /»>0.05). The mean Iq was 
higher in the MED than in the GMX (Mann-Whitney, 
J7 = 0.0295), but this difierence was only significant in June. 
Within the GMX, the Ig for females sampled in June was 
significantiy lower than for those sampled in the other months 
(Kruskal-WaUis,j7< 0.0125; Table 2). 

Histology 

Histological analysis of ovarian tissue from the GMX sampled 
from April to June showed no significant difiFerences in the gonad 
development between years. The number of inactive (IN) femEiles 
was less than 20% of our sample except in 2007 (28.0%). The 
proportion of active non-spawning (ANS) females was 20.0% 
(2007), 50.0% (2008) and 30.0% (2009). The proportion of active 
spawning (AS) females was consistent at about 50% for all years 
(Figure 2). Given the lack of annual variation, GMX samples were 
pooled for aU years for subsequent analyses. 

When ovary samples were arranged by month, consistent 
differences in the reproductive condition were observed with 
increasing maturation throughout the sampling period (Figure 3). 
In the GMX, the samples coUected in February and March did not 
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Figure 1 . Stages of oocyte development observed in Atlantic bluefin tuna sampled from the Gulf of Mexico and Mediterranean Sea. 

A: lipid stage follicles (LSF); B: LSF, and early/late vitellogenic stage follicles (VF); C, D: migratory-nucleus follicles (MNF); E: hydrated follicle (HF); F, G: 

post-ovulatory follicles (POF); H: alpha (aAF) and beta(/iAF) atretic follicles. All scale bars are 100 ^im. 

doi:10.1371/journal.pone.0098233.g001 



include any AS females. The proportion of AS females increased 
from April to June until a small decline was observed near the end 
of June. Near the end of the presumed spawning season in the 
MED (mid-June-mid-July), proportions of AS, ANS, and IN fish 
were similar to those observed in the GMX in May (Figure 3). 

Spawning frequency 

Spawning fraction in the GMX was estimated to be 0.45 
throughout the sampling period as 49 out of 108 mature females 
contained POFs in their ovaries. When this parameter was 
calculated by month, the proportion of females with POFs caught 
in April was lower than that in May and June. When the spawning 
frequency was calculated considering only AS females, its value 
increased significantly and remained similar among months. The 
spawning fraction in the MED samples was higher (0.60) than in 
the GMX (Table 3). 

Stereology 

No significant differences were found between months for the 
stereological parameters of atretic follicles (aAF and /^AF) in fish 



sampled in the GMX (Table 2). No significant differences were 
found for the numerical density of LSF (A'l/LSF) in the GMX 
despite the decrease observed throughout the reproductive season. 
Nevertheless, the total number of LSF (jVLSF) and the relative 
number of LSF (jVgLSF) estimated in GMX fish in April was 
significantiy higher than in June. In general, stereological counts of 
LSF quantified in the MED showed significantly lower values than 
those from the GMX sampled in April and May. Although the 
mean number of VF per mm^^ (NyYF) and the mean number of 
VF (jWF) remained unchanged in the GMX, the relative number 
of VF (.A(§VF) was significantly lower in June than in May. NyW 
and jWF estimated in the MED were similar to the GMX values 
in June; however, because eastern fish were smaller than western 
fish on average, the relative number of VF (.A^VF) was much 
higher for MED fish. The low number of females with MNF and 
HF is likely the cause for finding no significant differences for these 
stages throughout the sampling period in the GMX. The highest 
number of POF (TV^POF, jVPOF, and MgPOF) occurred at the 
beginning of the samphng period (April) and the lowest values 
were observed in June; however, no significant differences were 
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found either among months or between the spawning areas. The 
reahzed fecundity estimated as the mean absolute number of POFs 
( + SD) was 7.65-10''( + 6.7M0'') for fish from the GMX. This 
corresponds to relative batch fecundity (Ng^^) of 
28.14( + 26.90)POF g-'. In the MED, the estimated mean 
number of POF was 7.36-10''( + 6.7M0'^) and the relative 
fecundity 45. 56( + 33. 79)POF g"' (Table 2). 

Discussion 

Although not spatially and temporally exhaustive, this study 
represents the first attempt since the 1970s to accurately assess tiie 
spawning condition of ABFT sampled on the GMX spawning 
grounds and is the first histological and stereological comparison of 
the eastern and western spawning stocks. 

Histological analysis of gonad samples from the spawning 
grounds throughout the spawning season is essential for evaluating 
the reproductive condition and performance of ABFT. Addition- 
ally, systematic sampling across the extent of the spawning 
grounds allows the study of temporal variations in key reproduc- 
tive parameters, such as sex ratio, proportion of mature fish, 
spawning frequency, and spawning periodicity and fecundity. This 
information was lacking for the western stock resulting in large 
uncertainties for stock assessment and evaluation of productivity 
[8]. Since the implementation of the moratorium on directed 
fishing for ABFT in the GMX, federal fisheries observers have 
sampled ABFT caught as bycatch in the yeUowfin tuna and 
swordfish longline fisheries {Thunnus albacares and Xiphias gladius, 
respectively). Such restrictions prevent comprehensive size sam- 
pling of spawning ABFT in the GMX and hinder the determi- 
nation of the spatial and temporal extent of western spawning 
[31]. 

As a result of the bycatch sampling, previous studies have lacked 
small/ medium fish ( < 1 80 cm) leading to larger/ older size and age 
at maturity estimates [20,21]. Despite electronic tagging data that 
show presumed mature fish outside of known spawning areas 
[1,32-34], these studies assumed western bluefin only spawn in the 
GMX. A New England study suggested previous ABFT sampling 
did not accurately represent the spawning size range of the western 
population because it only included fish sampled by longhners on 
known spawning grounds rather than all size classes sampled 
throughout their range [14]. Gear type, size selectivity, and 
vertical distribution of tuna by size also influence the size of 
spawners sampled by commercial fishing fleets [12,35]. 

In this study, the smallest tuna sampled from the GMX, 1 72 cm 
and estimated age 7-8 y [36], had ripe ovaries with numerous 
recent POFs. This is consistent with an earlier study finding 
mature fish at 8 y [14] but not the recently proposed 12-16 y 
[20,21]. In order to fully understand the reproductive dynamics of 
the western spawning stock, the maturity ogive should be revised 
using comprehensive size samphng over larger temporal scales 
including histological examination of the ovaries and endocrine 
proffling. New studies utflizing endocrine hormones developed and 
calibrated in captive eastern ABFT [37] indicate that western 
ABFT become sexuaUy mature at less than 8 y [38]. 

Prior to the US moratorium, catch records indicated the 
presence of only giant ABFT (>180 cm) in the GMX [3]. As 
opposed to the current management paradigm of western ABFT 
maturing at an older age than the eastern stock, fish may exhibit 
size and temporal segregation on the spawning grounds as is seen 
with the eastern spawning stock [3,32,39] and in Pacific bluefin 
tuna [Thunnus orientalis) [40] . There is indirect evidence that smaUer 
fish may utflize alternative spawning locations, such as the 
Caribbean Sea, the Bahamas, or the Gulf Stream margins because 
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2007 2008 2009 



Figure 2. Maturity stages for female Atlantic bluefin tuna 
sampled in the Gulf of Mexico from 2007-2009. No significant 
difference was found between years for maturity stages, and tinus, all 
years were pooled for future analyses. IN = inactive, ANS = active non- 
spawning, AS = active spawning. 
doi:1 0.1 371 /journal.pone.0098233.g002 

ripe fish have been sampled there [9,15,41]. Given this evidence, it 
is possible that smaller bluefin spawn in alternative locations 
within the GMX as previous sampling has been concentrated in 
the north/central Gulf where US longline vessels operate. An 
ABFT life history model predicts that smaller /younger maturing 
fish should have shorter migration routes and spawn in areas closer 
to feeding areas than larger, older fish with higher energy reserves 
[42]. 

Electronic tagging results have consistendy shown annual 
migration patterns of giant ABFT not entering either known 
spawning ground before returning to northern feeding grounds 
[1,32-34]. It is possible that western ABFT spawn over a broader 
area of regions with oceanographic conditions appropriate for 
larval development than previously assumed [3,14,32,34,43]. 
Recent larval cruises found bluefin larvae outside the GMX 
[44], but spawning areas beyond the northern GMX await 
histological validation. Pop-up satellite tagged 2-5 year old ABFT 
did not enter the GMX or MED during presumed spawning 
periods (April-June) [45]. Nevertheless, some fish lingered in 
subtropical seas north of the Bahamas and in the southern mid- 
Atlantic Bight, areas visited by tagged adults [1,5,33,46]. 

Our histological examination of the GMX bluefin ovaries 
revealed differences in follicle maturation between months 
throughout the spawning period. While samples collected in 
February and March contained no active spawning (AS) individ- 
uals, 3 1 % of samples collected in April were AS individuals. As the 
spawning season progressed, the number of AS females increased 
and peaked in May (60%). The /g observed in June was 
significantiy lower than other months indicating a decrease in 
ovarian size, and thus, the impending end of the spawning season 
for this region. The statistical results of the stereological analysis 
are consistent with previous findings and with the progression of 
the spawning season[l 1]. As the spawning season progresses, LSFs 



GMX MED 




Feb Mar Apr May Jun Jun/Jul 



Figure 3. Maturity stages for female Atlantic bluefin tuna 
sampled in the Gulf of Mexico and Mediterranean Sea 
separated by month. No significant difference was found between 
years for maturity stages, and thus, all years were pooled for future 
analyses. IN = inactive, ANS = active non-spawning, AS = active spawn- 
ing. 

doi:10.1371/journal.pone.0098233.g003 

become less frequent indicating high levels of recruitment to VFs 
thereby compensating for losses caused by atresia or spawning. 
Similarly, the relative number of VFs {NgYF) was significandy 
higher in April and May than in June indicating a decrease in the 
recruitment of VFs as the end of the spawning season approaches. 

In spawning fish, atresia is a natural mechanism for regulating 
the number of eggs spawned. Alternatively, massive atresia can 
indicate a cessation of oocyte maturation and/or spawning 
activity [47]. The GMX samples showed relatively low and stable 
levels of aAF throughout the spawning season indicating ABFT 
found in the GMX are in favorable condition for oocyte 
maturation and spawning. While these fish appear to be actively 
spawning, tagged ABFT of presumed reproductively mature size 
observed outside the known spawning areas during the spawning 
period could be skipping spawning due to unfavorable body 
condition [14,34]. ABFT sampled on the New England and 
Canadian foraging grounds have had periods of reduced somatic 
condition [48,49] possibly accounting for an increased incidence of 
skipped spawning [50]. The incidence of skipped spawning in 
ABFT, however, is unknown, and modeling results show it is less 
likely to occur in larger, older fish in positive energy balance [42] . 
Giant ABFT sampled on western foraging grounds in the fall 
[14,48] and in the MED in the spring and early summer [51] have 
extensive perigonadal fat and somatic lipid stores, and thus seem 
unlikely candidates for skipped spawning. 

While ABFT have been observed on the western spawning 
grounds for months during the spawning season [34], individuals 
are believed to be actively spawning for only a few weeks 
[9,52,53]. These findings, albeit fishery-dependent, define the 
temporal borders of the reproductive events occurring in the 
north/ central part of the GMX indicating the spawning season 
runs from April to June with maximum spawning activity in May. 
However, ABFT begin entering the GMX in late November, and 
those arriving in winter experience warm water masses of >24°C 
in the lower GMX [34]. Reproductive sampling has been 
primarily conducted in the northern GMX and US territorial 
seas in late spring [9,52]. Spawning activity occurring earlier in 
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other areas of the GMX awaits confirmation by broader sampling, 
especially in Mexican territorial seas. 

The relatively low proportion of GMX females with POFs in 
their ovaries (<51%) contrasts with the high spawning frequency 
(60%) observed in the western MED [12,13]. The lower spawning 
frequency observed in the GMX could be the result of bias 
associated with utilizing the yeUowfm and swordfish fisheries as the 
only sampling method. As long as bluefin reproductive studies rely 
on bycatch in commercial fisheries, it is not possible to obtain an 
unbiased, accurate assessment of ABFT reproduction. Given these 
constraints, it is important to note the temporal and spatial aspects 
of the sampling as well as the fishing gear used for any bluefin 
maturity study. 

Stereological methods have often been used as an accurate tool 
for estimating fecundity in fishes, including eastern ABFT 
[11,22,29,54,55]. Realized fecundity can be estimated through 
stereological counts of POFs, whereas the number of MNF is an 
estimation of the potential fecundity [29,56]. In this study, the 
mean relative batch fecundity was calculated directly from 
stereological counts of POFs and showed a decrease as the season 
progressed. This is atypical for indeterminate spawners [57,58] but 
is likely due to the selective nature of sampling bluefin as bycatch 
in a longline fishery. One study suggested that monthly variation 
of fecundity may be masked by a decrease in the condition factor 
of fish appearing later on the spawning ground [59]. Although 
significant differences were not found between months, the highest 
value in the relative fecundity occurred early in the season (April), 
even though the number of AS females was still quite low. ABFT 
entering the GMX early might exhibit higher reproductive 
potential than those arriving later due to the good condition 
acquired on the foraging grounds. Otherwise, the lower spawning 
frequency observed on the western spawning grounds could be a 
consequence of migration distance [42] or decreased body 
condition observed on the western foraging grounds [48]. While 
being the first to arrive on the spawning grounds might provide 
increased resource availability for offspring, arriving in poor 
condition could decrease larval survival rates [51,60]. 

The fecundity of eastern spawning bluefin was estimated at 
59 eggs-g-' [12] and 48.22 eggs-g"' [13] for potential and 
realized fecundity, respectively. Our results show a lower realized 
fecundity for western spawners (28 eggs-g^') than eastern spawn- 
ers; however, this difference was not statistically significant. This is 
in agreement with previous work indicating realized fecundity is 
not proportional to body size [13]. Although the differences in 
fecundity and maturity schedules of both ABFT stocks can be 
masked by sampling bias and complex population dynamics, 
fecundity has been shown to vary within a given species as a result 
of different adaptations to environmental habitats [61]. The 
unexpected variation of fecundity shown among large ABFT could 
be related to the environmental adaptations and the balance 
between body condition and energetic expenses during migration 
E „ [42,51]. 

Eastern and western ABFT spawning sites seem to exhibit the 
same periodicity (three months), but spawning in the northern 
GMX occurs one month earlier than in the western Mediterra- 
nean spawning ground. This is likely due to specific oceanographic 
conditions and the early warmer temperatures observed in the 
GMX (data from CCA-UNAM). In this study, we have observed 
similar values in bluefin reproductive parameters showing that the 
spawning condition of Mediterranean spawners from mid-June to 
mid-July is comparable with the reproductive peak observed in the 
GMX in May. Depth and temperature associations of electron- 
ically tagged ABFT entering the GMX in winter also suggest that 
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it also serves as a foraging ground [34], as it does for other tunas 
and billfish species. 

A more holistic view of the population dynamics of ABFT 
requires life history characteristics, reproductive profiles, and 
spawning areas and periodicity to be well defined, especially since 
they undoubtedly will change with shifts in climate and ocean 
productivity [62] . The extent and quality of lipids acquired by 
tunas before they arrive in spawning areas will affect eggs and 
larvae, and therefore, overall stock biomass [42,v^l]. Future work 
should address energetic relationships between reproduction, 
migration, and early life history through modeUing, biological 
sampling, and the development of smart tags to detect actual 
spawning events. 
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